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1 Introduction: 

Given that it’s an exciting news topic due to the studies that took place not so long 
ago regarding the railway line that’s going to connect MARRAKECH, AGADIR and RABAT 
I thought that it would be interesting to research the optimization of the railway profile. 
wheel tread wear is a form of wheel damage that can seriously affect the performance of 
freight vehicles. A new numerical approach to optimizing wheel profiles can reduce 
circular wear on the LM wheel in the design cycle. This approach considers the influence 
of different line conditions and speed fluctuation on wheel wear, along with the 
performance of the wheel and the rail as the materials wear. In this approach, a nonlinear 
numerical optimization model for the wheel tread profile is built through a 
backpropagation (BP) neural network method. 


The multipoint Kik-Piotrowski (KP) contact mechanics model is applied to calculate the 
wheel/rail normal force, tangential creep force, the stick-slip area, and the size and shape 
of the contact patch. The optimal profile is obtained through the genetic algorithm (GA) 
method. In order to better reflect the random characteristics of wheel/rail matching and 
interval uncertainty, a random sampling technique is used to generate a random data 
sample at typical operating speeds. 


This report proposes a rail profile optimization method that takes account of wear 
rate within design cycle so as to minimize rail wear at the curve in heavy haul railway and 
extend the service life of rail. Taking rail wear rate as the object function, the vertical 
coordinate of rail profile at range optimization as independent variable, and the geometric 
characteristics and grinding depth of rail profile as constraint conditions, the support 
vector machine regression theory was used to fit the nonlinear relationship between rail 
profile and its wear rate. Then, the profile optimization model was built. Based on the 
optimization principle of genetic algorithm, the profile optimization model was solved to 
achieve the optimal rail profile. A multibody dynamics model was used to check the 
dynamic performance of carriage running on optimal rail profile. The result showed that 
the average relative error of support vector machine regression model remained less than 
10% after a number of training processes. The dynamic performance of carriage running 
on optimized rail profile met the requirements on safety index and stability. The wear rate 
of optimized profile was lower than that of standard profile by 5.8%; the allowable 
carrying gross weight increased by 12.7%. 
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2 Genetic algorithm: 

GAs are based on the principle of genetic evolution in which only the best or “fittest” 
individuals of a population survive to reproduce and create the next generation. During the 
process, some genes of the new-born individuals are randomly mutated. This operation 
introduces new changes that could make an individual stronger to survive in its environment, 
and therefore, more likely to reproduce and pass those advantages to the next generation. In 
the mathematical analogy a population consist of several design vectors (individuals), each 
design vector is composed of a set of design variables (chromosomes) that can be modified 
within a specific domain or search space during the algorithm execution. In the case of single- 
objective optimization problems, a single objective function that assigns a fitness value to the 
design vectors is defined. The natural processes of selection, reproduction and mutation are 
simulated through numerical operations. At every iteration (or generation), the best design 
vectors of the population based on their fitness value are selected, combined, and mutated 
to create a new stronger population. The algorithm usually stops after a certain number of 
generations is completed, if no improvement is observed after a specific number of iterations 
or if another convergence criterion is satisfied. 


3 Positioning of the problem: 

The rail profile is considered invariable and free from the impact of wear-induced 
change in rail profile during wheel-rail contact in respect of rail profile optimization. The 
rail optimization result is usually the current simulation comparison result for 
demonstrating new wheel-rail profile matching effect without considering the fact that 
rail profile is subjected to wear all the time during operation, that the rail profile changes 
all the time due to wear, and that the worn wheel-rail profile contact relationship changes 
therewith. From this point of view, the current wheel-rail contact relationship in optimum 
state is not in the same position to assure the optimal matching effect after the operation 
of train for a period of time. Hence, this paper proposes a rail profile optimization method 
that takes account of wear rate within design cycle. The maximum wear depth, that is, 
1mm, was taken as design cycle. This method updates rail profile through wear loss 
calculation based on random irregularity on rail with optimized profile as the initial 
profile, subjects the updated rail profile to another dynamics simulation, calculates the 
wear loss of updated rail profile, and then works out the wear rate. Since the final profile 
is immediately associated with initial design profile and the contact geometry state during 
the operation of train, the paper established the nonlinear relationship between wear rate 
and rail profile geometric parameter within design cycle through nonlinear fitting, built 
the optimization model, sought to solve the model with optimization algorithm, and 
sought an optimized rail profile with the minimum wear rate. The optimization process is 
shown in Figure 1. 
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Figure 1: Flowchart of rail profile optimal design. 


4 Mathematical model for rail profile optimization: 

This paper studied the regular pattern of rail wear at small-radius curved section 
of heavy haul railway and established a numerical model for profile optimization with the 
high rail profile at curve radius of 600 m along a chosen heavy haul railway. Assuming that 
the coordinate of rail head cross-sectional point is (xi,yi) and the object function is wear 
rate F(x,y), then the rail profile could be expressed as spline interpolation 
function f(x,y) of rail head coordinate point. The mathematical model between rail wear 
rate and profile is nonlinear mapping relationship; let 
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Subject the mathematical model to optimal solution to get the rail profile with the 
minimum wear rate. 


4 1 Determination of Independent Variable : 

According to research findings, premature renewal rail is necessary for CHN 
75 kg/m rails in general use for heavy haul railways in China due to side wear when the 
carrying gross weight reaches approx. 350 Mt (2013) [20]. Figure 2 [4] shows the wear 
evolution of high rail profile at curved section of a heavy haul railway with curved section 
radius of 600 m designed with standard CHN 75 kg/m rails. 
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Figure 2: Diagram of rail head optimization range (unit: mm) 


As shown in Figure 2, the seriously worn points of high rail are principally located 
between rail center line and gauge side. Hence, the points within rail profile optimization 
interval correspond to horizontal coordinate -10~36 mm. 12 discrete points are chosen 
within optimization interval as optimization points. In view of the fact that rail wear 
principally changes the vertical coordinate of rail profile point, the independent variable 
in mathematical model for rail profile optimization should be: 


= V2 esse y12)T 
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4.2 . Determination of Dependent Variable : 

Rail wear is a process of accumulation during operation of carriages. The wear loss variation 
mechanism of different rail profiles is not constant within overall interval for different 
segments of a time interval. Since the grinding for rail profile optimization is performed within 
a certain cycle, the rail wear loss shall be changed at determination of object variable so as to 
avoid local optimum. Expressing the service efficiency of rail as rail wear rate within section is 
a better way to reflect the wear accumulation process. Take T as the time interval of design 
cycle, and divide T into n intervals, Mi is the rail wear within the ith interval, Mi represents 
the number of carriages passing within the ith interval, and then the wear rate could be 
expressed as: 


r= yn i=1 mi/ yn i=1 Mi. 


Take the unit length of rail at maximum wear depth as the object of study, and the 
accumulative wear loss would be: 


> mi = [S0 (yn 1, yn2.,...,yn 12)-ST (yO1,y02.,...,y012)] - L, 


where the rail length lis equal to 1, that is, the unit length. SO and ST, respectively, represent 
the area between initial rail profile and horizontal coordinate axis and the area between rail 
profile and horizontal coordinate axis after time T. 


4.3 Determination of Constraint Condition : 
According to wheel-rail contact relationship and the pertinent literature (2009), the 
determination of rail profile must meet convex curve condition; that is to say, the slope of 
segment between consecutive points decreases with the increase of horizontal coordinate. 
The following constraints shall be imposed for vertical coordinate of the 12 optimization 
points: 


(yi - yi-1)/ (xi - xi-1) - (yi+1 - yi-1)/ (xi+1 - xi-1)>0(i=1,...,12). 


Since rail profile is optimized based on original standard profile, the optimized profile is below 
the original profile, while the optimization of rail profile shall meet the requirements on max. 


The mathematical model established for rail profile optimization based on analysis of each 
element of single-object optimization model is as follows: 


TT  : 
2020/2021 7 








au 


LAL yl pull 44 dl 





Caw digaU Adored! dv roll 


Ecole Mohammadia d'Ingénieurs projet techniques d’optimisation te ir ns 
min r (y1, y2,...,y12) = Yn i=1 mi /¥n i=1 Mi 


s.t. ai < yi < bi 
0<Ayi<3 (i=1,...,12) 


D (y1, y2,...,y12) < [d], 


5 PSO-SVM Regression Fitting Rail Wear Rate Model: 


There is a complex nonlinear relationship between rail wear rate and rail profile. The 
dependent variable and independent variable are subjected to nonlinear regression through 
given sample data based on principle of statistics. A mathematical model with generalization 
capability is attained through sample training for machine learning. Machine learning is 
designed to handle any complex nonlinear relationship. It converts input and output data by 
seeking the “black box” of hidden layer, but the complex model obtained through training 
cannot be expressed with an explicit relation. This paper identifies the nonlinear relationship 
between rail wear rate and rail profile in this way. 


5.1 The Sample of Rail Wear Rate: 
The sample for data fitting (the wear rate of different rail profiles) can be obtained by field 
measurements and multibody dynamics simulations. Due to the number of rail wear 
observations, they are unable to meet the sampling requirements. Therefore, this paper used 
a multibody dynamics simulation to obtain the rail wear rate. 


The main purpose of sampling is to analyze the selected samples and to deduce the 
characteristics of all the samples. In order to obtain more comprehensive rail wear rate data, 
this paper takes a sampling method to obtain the samples needed for data fitting. According 
to the characteristics of rail wear and grinding maintenance, the rail profiles were sampled. 
As the wear of the rail profile is a continuous change, in order to make the sample profile 
representation and uniformity and other characteristics, the D-optimal experimental design 
sample was selected. According to the composition of the curve to be smooth and smooth by 
the optimization points of rail profile and the geometrical requirements of the rails in the 
actual railway condition, the D-optimal experimental design was used to sample 101 sets of 
rail profiles. Based on the vehicle-track coupling dynamics model and the Archard wear model, 
the wear of the sample rail profiles was simulated and the wear rate of the rail profiles was 
obtained. 


5.2 SVM Nonlinear Regression : 
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Support vector machine, also known as SVM, is a learning method where computer is used for 
statistical analysis of data. The rationale of SVM is that the sample relationship identified 
through summarization and generalization of given sample point data is generalized. The basic 
means of SVM regression analysis for rail wear rate model is linearization that converts the 
nonlinear problem between profile and wear rate into the combination of linear problems in 
higher dimensional space. In consideration of the rationale of SVM, assuming that y = (y1, 
y2,...,y12)" is input sample vector and r is the corresponding output value, then: 


r=f(y)+V 
where V is the offset parameter to be determined. 


Seek the optimal regression surface analysis rationale based on SVM, define two nonnegative 
slack variables (i.e., č and 7), constant penalty factor C, and insensitive error parameter £, and 
realize the conversion into a quadratic convex programming problem through regression of 
sample data. 


min 1/2 ||w||? + C m > i=0 (či + ni) 
s.t. ri (w - yi)-bse+éi 


(w: yi)+b-rise+ni (i=1,2,...,m). 











Use the Lagrange multiplier method and introduce operators ai and ali; then the optimal solution 
to the rail profile optimization problem in higher dimensional space is: 


fly)= (w - y) +b 
= ÿ support vector (ai - af i) K(y, yi) + b, 





where K(y, yi) = (p(y) : g{vi)) 


5.3 SVM Regression Parameter Optimized with PSO: 

penalty factorCand kernel function radius g significantly affect the regression prediction 
accuracy of model during SVM-based model training. Under normal conditions, C is 
monotonously negatively correlated with fitting error within a certain range. Where the same 
value is chosen, prediction error nonmonotonically decreases with the increase of penalty 
factor, but the increase of its value may result in overfitting. Where g value is small, the kernel 
function curve changes slowly, in which case overfitting may arise during model fitting and 
affect the prediction; underfitting may arise when g value is small. Hence, particle swarm 
optimization (PSO for short) is used for optimization and to identify appropriate C and g values 
during regression of rail profile optimization model. Mean square deviation is taken as the 
fitness function of model when SVM parameters C and g are optimized with particle swarm 
optimization. The obtained relationship between best particle means square deviation and 
the number of iterations is shown in Figure 3. 
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FIGURE 3: Swarm mean curve of particle swarm optimization. 


As shown in Figures 3 and 4, swarm mean square deviation gradually decreases with the 
increase in the number of iterations within a certain range. When the number of particle 
iterations is less than 200, the swarm mean square deviation curve changes significantly. 
When the number of iterations reaches 300, the training swarm mean square deviation curve 
gets into a convergence state and converges at 0.013. The best particle mean square deviation 
curve substantially tends to be stable after 50 iterations, and the convergence value is 0.012 
795. The corresponding optimized kernel function radius g = 0.006 8, and the penalty factor 
C= 101.866 4 
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FIGURE 4: Optimal value curve of particle swarm optimization. 
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5.4 Rail Wear Rate Model Analysis: 

After the optimal kernel function radius g and penalty factor C are obtained through particle 
swarm optimization, K cross validation is used to 0.0126 0.0128 0.0130 0.0132 0.0134 0.0136 
Best particle mean square deviation 0 50 100 150 200 250 300 Number of iterations Figure 4: 
Optimal value curve of particle swarm optimization. Real SVM 0 2 4 6 8 10 12 14 16 18 20 Rail 
profile test set number 0 2 4 6 8 10 High rail wear rate (1 x 10-5 mm2/t) Figure 5: Rail profile 
test set prediction result of wear rate. subject a certain part of the 101 selected groups of rail 
profile wear data to random sampling of rail profile. The rail samples taken are divided into 
training set and test set. Training set is used for modeling during regression, while test set is 
used to determine the regression accuracy of model. The stability of built model is analyzed 
and compared through a number of cross combinations. In this paper, 20 groups of data 
constitute the test set for model; whole others constitute the training set for the regression 
of rail profile optimization model. Figures 5 and 6 show the prediction result of test sample by 
PSO-SVM, and Table 2 shows the prediction data of test set by PSO-SVM rail profile 
optimization model. 
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FIGURE 6: Rail test sample prediction wear rate error. FIGURE 5: Rail profile test set prediction result of wear rate 





As shown in Figures 5 and 6, the predicted value of test sample is extremely close to its actual 
value, and the variation tendency of predicted value substantially keeps up with that of actual 
value. According to the prediction data in Table 2, the overall relative error of prediction for 
test sample is less than 20%. The prediction errors for Groups 40 and 22 are beyond 10%, 
while the prediction errors for all other 18 groups of samples are around 5%. The mean 
relative error of prediction for the 20 groups of test samples is 4.3%, which means that PSO- 
SVM regression model for rail wear rate has excellent prediction performance. 


6 Numerical Solution to Rail Profile Optimization Model: 
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Rail profile optimization involves a number of factors. Rail wear is associated with its material, 
the acting force between wheel and rail and other wheel-rail system factors, and is closely 
associated with ambient environment and subsequent maintenance and grinding. Upon the 
establishment of rail wear rate model, the optimal solution shall be found for rail wear rate 
model so as to identify the rail profile with minimum wear loss at curved section within design 
cycle with the purpose of reducing rail wear at curved section and extending the service life 
of rail as stated in this report. 


With the vertical coordinate of point at rail profile as the independent variable and the wear 
rate as the dependent variable, this paper identifies the randomness, complexity, multiple 
constraints, and so forth of nonlinear rail profile optimization model for curved section by 
subjecting both variables to fitting through PSO-SVM. Genetic algorithm is used to seek 
numerical solution of this model based on above noted features of the model. 


6.1 Parameter Setup of Genetic Algorithm: 
The following issues shall be taken into account when using genetic algorithm to seek to solve 
rail profile optimization model. 


(1) Swarm and Swarm Size. The independent variable y in rail profile optimization model 
represents gene, so the individual rail profile could be expressed as yi = (y1i, y2i,...,y12,i) T, 
where the aggregate of rail profile yi represents the swarm stated in this paper. Generally, the 
size of swarm is directly proportional to the all-sidedness of genetic algorithm solution, but 
larger size of swarm may consume longer time of iterative computation. The swarm size is 
defined as 200 based on requirement of this report. 


(2) Selection Operator. According to the “survival of the fittest” principle in genetic algorithm, 
the operators with more favorable fitness values are more likely to be chosen; otherwise, the 
individuals will be knocked out. Genetic operators in common use include proportional 
selection, the best individual preservation strategy, and deterministic sampling, and this paper 
takes proportional selection as selection operator. 


(3) Crossover Operator pc. Where genetic algorithm is concerned, “crossover” means genetic 
recombination; that is to say, two individuals exchange gene with each other in a certain way 
to generate a new individual. This paper performs two-point crossover for chromosome in the 
following procedure. Assume that y1 = (y1 1 , y1 2 , y1 3 ,...,y1 11, y1 12) T and y2 = (y2 1, 
y22,y23.,...,y2 11, y2 12) T; if gene exchanges locations r1 = 3 and r2 = 11, the chromosome 
after two-point crossover should be yR 1 = (y1 1, y1 2 , y2 3 ,...,y2 11, y1 12) T and yB 2 = (y2 
1,y22,y13.,...,y1 11, y2 12) T. 


(4) Mutation Operator pm. Mutation offers the genes omitted during computation, enables 
the global search of rail profile, and helps to avoid precocity. Single-point mutation is 
employed in rail wear profile optimization model. The operating procedure is as follows. For 


chromosome y = (1.1, 1.5, 2.2, 2.3, ..., 13)T, the mutation position is 2, and then yA = (1.1, 
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1.8, 2.2, 2.3,..., 13)T is obtained after random mutation. The process of seeking optimal rail 





profile through genetic algorithm is shown in Figure 7. 


Rail profile vertical coordinate coding 


Vertical coordinate combination 


Vertical coordinate mean square error 


= Yes 


Perform selection 
Rail profile crossover Rail profile decoding 


Single-point mutation Output vertical coordinate 


Generation of subgeneration rail profile 


FIGURE 7: Computation flow chart of genetic algorithm. 









6.2 Rail Profile Optimization Result : 
Figure 8 shows the optimized profile Opt that meets constraint condition and is computed 
by seeking to solve PSO-SVM rail profile optimization model through genetic algorithm. 
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FIGURE 8: Comparison between optimized profile and standard 
profile. 


Figure 9 shows the variation of vertical coordinates of the optimized profile and standard 
profile. Thus, it can be seen that the vertical coordinate difference at optimization point 
increases with the increase of horizontal coordinate. When the horizontal coordinate of 
optimization point is (-10, 20), the vertical coordinate changes slowly (normally no more than 
0.5 mm) at points in front and behind this point; when the optimization point is (20, 36), its 
vertical coordinate changes significantly, and the vertical coordinate difference increases with 
the increase of horizontal coordinate before and after optimization. The rail profile changes 
significantly when horizontal coordinate is located between 20 mm and 36 mm, and the 
maximum depth of Opt and standard profile optimization point variation is 1.900 4 mm. This 
is principally because the overall linearity of rail head exhibits a convex curve. Since the most 
frequent rail wear at curved section is side wear, the optimization depth at rail side shall be 
greater than at top face in order to reduce the contact between rail side and wheel flange and 
concentrate the wheel-rail contact at top of rail head to the greatest extent so as to minimize 
rail wear. 
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FIGURE 9: Local comparison of optimization interval (unit mm). 
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7 Analysis of Performance of Optimized Rail Profile: 


7.1 Wear Analysis : 
Genetic algorithm solution seeking and simulation calculation show that the optimized rail 
profile offers lower wear rate and improves the allowable gross rail load on axle as compared 
with standard profile under the same conditions within an established wear cycle, having 
reached the research objective of this report. 


In view of the progressive accumulation of rail wear, 1 mm of the profile wear depth was taken 
as rail wear analysis cycle; that is to say, the termination condition for simulation of rail profile 
wear was that the maximum wear depth reached 1 mm. In this range, the variation of the 
optimized profile wear determined through rail wear model simulation calculation, shown in 
Figure 10. Opt-0 is the optimized rail profile, and Opt-10 is the rail profile when wear depth 
reaches 1 mm. Figure 10 shows the accumulative variation of rail optimized profile wear 
depth. 











Vertical coordinate of rail profile (mm) 





-40 -30 -20 -10 0 10 20 30 40 


Horizontal coordinate of rail profile (mm) 


--- Opt-10 
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FIGURE 10: Profiles before and after wear of the optimized profile. 


It is observed from Figures 10 and 11 that most wear of high rail before and after optimization 
is located at horizontal coordinate 10-30 mm. Wear simulation calculation showed that the 
allowable number of passes at the optimized profile Opt is 707,605 vehicles, while the 
allowable gross load is 84.913 Mt, 12.7% up as compared with standard profile. The wear rate 
simulation value of the optimized profile Opt is 1.715 62 x 10-5 mm2 /t, 5.8% down as 
compared with standard profile. 
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FIGURE 11: Accumulation of the optimized profile wear. 


The wear rate of the optimized profile Opt predicted by genetic algorithm (PSO-SVM 
Regression) is 1.617 02 x 10° mm2 /t. Simulations comparison analysis showed that the 
relative error between rail wear prediction value and simulation value is 6%. To sum up, the 
result of optimization of rail profile with mathematical model of rail wear presented in this 
paper is consistent with optimization objective. 


7.2 Dynamic Performance Check for Optimized Profile: 

Carriage safety and stability constitute the prior condition for all maintenance and repair of 
rails. Dynamics analysis was performed on the motion state of a freight train with 30 t axle 
load passing through a curved section of 600 m in radius at 72 km/h. The safety and stability 
indicators like dynamic derailment coefficient, wheel-rail lateral force, wheel load reduction 
rate, and vibration acceleration of a train passing through a curved section were calculated 
with optimized profile as object of study under the same simulation conditions. Table 3 shows 
the maximum values of safety index for each wheelset in the case of optimized profile. 


TABLE 3: Maximum values of dynamics index of train in the case of optimized profile. 





Wheelset 1 Wheelset 2 Wheelset 3 Wheelset 4 
Safety index Standard values 
IL IR 21 2R 3L 3R 4L 4R 
Derailment coefficient 0.524 6 0,519 6 0.4015 0.403 3 0.564 8 0.375 0 0,337 2 0,472 0 0.8 
Wheel-rail lateral force/kN 93.169 59.631 60.741 68.655 70.991 41.896 42.635 52.804 180 


Wheel load reduction rate 0.637 0,541 0,489 0.58 0.8 








The body vibration acceleration curve of a train along the established route in the case of 
optimized profile is shown in Figure 12. 
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FIGURE 12; Vibration acceleration curve in the case of optimized profile 


The calculation values in Table 3 show that, according to the provisions of Dynamic 
Performance Test Appraisal Method and Evaluation Criteria for Rolling Stock (TB/T2360- 1993) 
on above-noted indicators [26, 27], the train complies with safety index during operation on 
optimized profile. According to the provisions of GB 5599-85 on vibration acceleration of 
freight trains that the vertical vibration acceleration shall not be greater than 0.7g (6.86 m/s?, 
g is the gravitational acceleration) while the lateral acceleration shall not be greater than 0.5g 
(4.90 m/s?), Figure 12 shows that the stability index of train traveling on optimized profile is 
satisfactory. 


8 Conclusion: 

This report established a mathematical model for rail profile optimization with purpose of 
reducing rail wear rate. During the mathematical modeling, the wear characteristics of curved 
portion of rail were taken into consideration while the 12 discrete points selected in 
optimization interval were used for spline interpolation to indirectly express the geometric 
profile of rail. To simplify the modeling process, unit-length rail was taken as the subject 
investigated, while the rail wear volume was translated into wear area, which could be 
calculated through the difference between profiles before and after wear. In such a manner, 
the independent variable and dependent variable of mathematical model were, respectively, 
numeralized. The following conclusion was reached: 


(1) Support vector machine was used for regression fitting of rail wear model based on sample 
rail profile simulation data. To improve the generalization capability of regression model and 
reduce the model prediction error, particle swarm optimization was proposed to optimize the 
kernel function radius g and penalty coefficient C in support vector machine regression. Rail 
wear rate model was built after sample training with PSO-SVM. The prediction of test set 
identified that the average relative error of test samples in this paper was 4.3%. After a 
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average relative error of test set was kept below 10%. 


(2) Genetic algorithm was used to solve rail profile optimal model to get the optimal rail profile 
based on optimality principle. According to simulation calculation, the wear rate of optimized 
profile Opt decreased by 5.8%, while the allowable gross load increased by 12.7% as compared 
with standard profile. 


(3) The analysis of train operation safety and stability on the optimized rail profile Opt 
performed with derailment coefficient, wheel-rail lateral force, wheel load reduction rate, and 
vibration acceleration as indicators showed that the optimized rail profile Opt met normal 
travel conditions of train. 


In consideration of all characteristics mentioned above, this mathematical model is available 
to optimize the design of curved portion profile of rail. 


This report takes into account the profile variation during service of rail and leaves out of 
consideration the diversity of wheel profile and the wear-induced change in wheel profile, 
which constitute the shortcoming of this study. To make up for such shortcoming, it is 
advisable to carry out in-depth research using the techniques similar to that stated in this 
paper, so as to obtain optimized rail profile with wheel-rail wear within the interval taken into 
account, thereby getting closer to practical situation. 


Comments: 


As for the implementation of this optimization method in Morocco, it is already in the process 
and engineers are adapting to these new methods since they collaborate with foreign 
companies that already benefited from them. The cost some times could be high, however, it 
turns out to be worth it in the long run due to the benefits and long-term upsides that these 
optimization methods provide. 
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